Spin Tunneling in Magnetic Molecules: Quasisingular 
Perturbations and Discontinuous SU(2) Instantons 

Ersin Kegecioglu and Anupam GargQ 

Department of Physics and Astronomy, 
Northwestern University, Evanston, Illinois 60208 
-- - (Dated: February 1, 2008) 

O ! Abstract 
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Spin coherent state path integrals with discontinuous semiclassical paths are investigated with 
special reference to a realistic model for the magnetic degrees of freedom in the Fes molecular solid. 
It is shown that such paths are essential to a proper understanding of the phenomenon of quenched 
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spin tunneling in these molecules. In the Fes problem, such paths are shown to arise as soon as a 
fourth order anisotropy term in the energy is turned on, making this term a singular perturbation 
from the semiclassical point of view. The instanton approximation is shown to quantitatively 

cn : 

explain the magnetic field dependence of the tunnel splitting, as well as agree with general rules 
for the number of quenching points allowed for a given value of spin. An accurate approximate 
formula for the spacing between quenching points is derived. 
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I. INTRODUCTION 



A. Overview 

Among the several advances in the study of large-spin molecular magnets in the last few 
years , perhaps the most remarkable is the observation of quenching of spin tunneling in the 
molecule [Fe 8 2 (OH)2(tacn) 6 ] 8+ (or just Fe 8 for brevity) by Wernsdorfer and Sessoli (WS) |2], 
|3| . The spin of this molecule can tunnel between classically degenerate orientations, but for 
certain values of the applied magnetic field, it is found that the tunneling is eliminated, or 
quenched. Between quenching points, the tunnel splitting oscillates as a function of applied 
field. The simplest model that is capable of describing this effect has been studied by many 
workers now j|, ||, |[ [7|, || || |10|, [Tl| , and it is well understood qualitatively in terms of an 
interference between tunneling Feynman paths. Such interference can also arise for massive 
particles in more than one dimension |12|. It arises in the case of a single spin because the 



kinetic term in the path integral has a geometrical or Berry phase structure j|, |l^, |l4jj . 

There is, however, yet another aspect of the Fe 8 experiments that has not been widely 
appreciated, because of which the physical explanation for the tunnel splitting oscillations 
is significantly incomplete. This is because if we view the tunneling in terms of just the in- 
terfering paths, the geometrical interpretation of the relative phase between them inevitably 
predicts a certain number — ten — of quenching points. Experimentally, and by numerical 
diagonalization of the appropriate model Hamiltonian, only four quenching points are seen. 
We will show in this paper that the resolution of this apparent paradox requires that we 
include paths that have discontinuities at the end points. As a general point, the necessity 
of including such paths in the spin-coherent-state path integral has in fact been known for 



some time [15], |16|| , but the significance of this idea has not progressed much beyond the 
purely formal level, and the one concrete problem where discontinuous paths are known to 
be needed — Larmor precession in a constant magnetic field — is so well understood to 
begin with, that few people have pursued this line of inquiry further. The Fe 8 problem has 
provided us with a strong motivation for doing just that, thus enhancing our understand- 
ing of the spin-coherent-state path integral. In turn, we can now see that discontinuous 
paths exist in all other coherent-state path integrals too, and this may prove to be of use in 
situations where such path integrals are the natural analytic tool. 
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There are, therefore, two points to this paper. First, as already stated, is to advance 
our understanding of the spin-coherent-state path integral. Second, we obtain a quantitative 
analytical explanation for additional aspects of WS's data that are not captured by the 
simple model, and that have up till now been only understood numerically. This is desirable 
as the oscillations seen in Fe§ are, to date, the clearest evidence for spin tunneling of a spin 
of such large magnitude. 

This may be a good point to ask why one should bother with an analytical explanation 



at all. After all, the model Hamiltonian (1A) given below can be numerically diagonalized 
with minimal effort. A numerical diagonalization by itself gives no insight into the results, 
however, and if it were the only way we had of solving the problem, the oscillations in the 
splitting would be a complete mystery. If one believes that the mental picture of interfering 
tunneling paths is useful, then it is surely important to know if and when and how this 
picture breaks down. Secondly, the numerical approach cannot by itself explain the scale of 
the tunnel splitting. 

A shorter paper with our results has been published previously [[17]]. In the present paper, 
in addition to providing the details of the work, we explain how to find the end points of 
the discontinuous paths or boundary jump instantons, as we refer to them. We also discuss 
the general model for parameter values other than those relevant to Fes in order to further 
test the scenario for the elimination of quenching points. Finally, we show that the number 
of quenching points that remain must be even or odd depending on the value of the spin, J, 
and how this comes about from the instanton analysis. 



B. Basic Facts about Fe§ 

The main facts about Fe§ are as follows. The cluster of eight Fe 3+ ions in one molecule 
has an approximate D 2 symmetry. Competing exchange interactions among the ionic spins 
in one molecule lead to a ground state with a total spin of J = 10. The compound forms 
a well ordered crystalline solid, in which the magnetic ions in one molecule are kept well 
separated from those in another molecule by large organic ligands. The magnetic degrees of 
freedom may thus be treated as a set of non-interacting spins, each with J = 10. To a first 
approximation, the anisotropy of each molecule may be described by the Hamiltonian 

n = -k 2 j z 2 + (h - k 2 )j x 2 - gf i B j ■ h. (1.1) 
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TABLE I: Important parameter combinations 



Quantity Formula 

H c 2k 1 J/gn B 

A k 2 /h 

A 2 CJ 2 /k x 

h H/H c 

C 4A 2 /i 2 

Here, g is the gyromagnetic ratio, hb is the Bohr magneton, and k\ > k 2 > are anisotropy 
parameters. Equation ( |1.1[ ) is the simplest Hamiltonian that can describe the quenching 
effect. Viewed as a classical Hamiltonian, it has two degenerate ground states at the points 
J = ±Jz in the absence of an applied field. These states are separated by a barrier k 2 J 2 
along the y-axis, and k\J 2 along the x-axis. Quantum mechanics admixes these states via 
tunneling [ 18| . When the field is non-zero, the classical minima are no longer along ±z, but 



they are still degenerate if H _L z, the easy axis, and one can have tunneling between the 
corresponding quantum mechanical states. More generally, one can have tunneling between 
excited states in the two wells. In this paper, however, we will only consider H||x, and 
tunneling between the lowest two states. In this case, the quenching points are perfectly 
regularly spaced, and are located at 

H x = (J - n - ±)AH(°\ n = 0,l,...,2J-l, (1.2) 

where the spacing or period is given by 



Atff = ^j^H c , (1.3) 

where A = k 2 /k\, and H c = 2k\ J/g^B- (These, and other important parameter combinations 
are tabulated in Table |.) 

We note in passing that although this and similar results for other quenching points were 
initially obtained by various semiclassical approximations 0, [5], 0], they are in fact exact (19 
for the model (II . II) . 



WS's observations differ from the simple model predictions in that they see only four 
quenching points on the positive H x axis, where the simple model would say ten, and the 
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FIG. 1: Tunnel splitting A between the ground level pair for the Hamiltonian flLJ ) with C = 
(solid line) and C = 29 fiK (dotted line). Plotted is log 10 |A|, with A measured in Kelvin. 

spacing between these points is more than 50% greater than the formula ( |1.3| ) would give. 
(See Fig. [I].) At the same time, this spacing is quite uniform. They surmise that the differ- 
ences are due to the presence of higher order anisotropy terms in the Hamiltonian consistent 
with the symmetry of the molecule, and by trial and error and numerical diagonalization of 
a 21 x 21 matrix, they discover that their data can be very well fit by the following model 
Hamiltonian [EOH 



H = -k x Jl + (ki - k 2 )J 2 £ - C[J\ + Jt\ - gVisJ • H. 



;i-4) 



The parameters g, k\ and k 2 are known through a variety of experimental evidence [21, 
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We will use the values used by WS: g ~ 2, k x ~ 0.338 K, k 2 ~ 0.246 K, and 
C = 29 fiK. It turns out to be important that one needs C > 0. 

This, then, describes the problem that we wish to solve. One point to note is that the 
dimensionless strength of the fourth order anisotropy C is 



A 2 = CJ 2 /h. 



1.5) 



For the Fes parameter values, A2 = 8.580 x 10~ 3 . That such a small term should have such 
a large effect (elimination of six quenching points, 50% increase in period), suggests that it 
is a singular perturbation. We shall see that this is indeed true when we analyze its effect 
on the semiclassical paths, although the perturbation is perfectly well behaved from the 
quantum mechanical operator point of view. 
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C. Plan of paper 



The plan of the paper is as follows. In Sec. || we briefly review the formalism of SU(2) 
instantons with special emphasis on the jump instantons. To this end, we first discuss 
(Sec. |il A| ) stereographic coordinates for the unit sphere, the corresponding parametrization 
of spin coherent states, and how paths on the complexified unit sphere are represented in 
terms of stereographic as well as the usual spherical polar coordinates. We then discuss 
(Sec. 11 .B) the action functional for spin, and show that unless one includes a "boundary 



term" that depends explicitly on the boundary values of the spin path, there is no sensibly 
formulatable least action principle, nor is there a solution to the Euler-Lagrange equations 
in general. Paths with discontinuities arise naturally once a boundary term is included. We 
also discuss Klauder's formulation of these points in terms of an extra infinitesimal pseudo- 
inertial term in the kinetic energy and the concomitant production of narrow boundary 
layers in the spin paths. In the limit where the extra kinetic energy term vanishes, the 
width of these layers also vanishes, and they then look like discontinuities in the paths. 

In Sec. |II C[ we discuss the general structure of the tunnel splitting for the model ( |1.4|) 
when both types of instantons — with and without jumps — are allowed. We shall see that 
an instanton of the latter type dominates as the field gets large, and since this instanton 
has no interfering partner, we see why the splitting ceases to oscillate after fewer than ten 
quenching points have been realized. 

In Sec. [TTI| we present the explicit analysis for our model for Fes, Eq. ( |1.4| ). In Sec. [Ill A 



we analyse the interfering instantons. We identify ( = AX 2 h 2 (with h = H/H c ) as the 
appropriate small parameter, and explain why the H x spacing of the quenching points is so 
regular. We obtain a concrete formula for this spacing, namely, 

where I is an integral expression [see Eq. ([3.271) 1 depending on the parameters A = k 2 /ki, 
A2 = CJ 2 /k\ (see also Table |I[). With the parameter values relevant to Fes, we obtain 
/ = 3.88, giving AH X = 0.409 T. WS quote a period of 0.41 T. Compared to the model Ho, 
the period is increased by a factor 7r(l — A) -1 / 2 // = 1.552. We relate I to a complete elliptic 
integral, and give various formulas for it in the Appendix. 



In Sec. I([B we turn to the jump instantons. We first formulate the equations for de- 



termining the endpoints of these instantons. In general these equations would have to be 



solved numerically, but for the Fes model, it turns out that an analytical solution can be 
found. With these endpoints in hand, we can find the complete instantons and their actions 
(this must still be done numerically). We use our computed actions for all the instantons 
to find the splitting as a function of magnetic field, and compare with the numerical diago- 
nalization. This exercise is carried out for the value of A2 applicable to Fe§, as well as some 
others. 



In Sec. [TV] we address the issue of the number of quenching points in light of general 
quantum mechanical theorems. It is not an accident that WS see an even number (four) 
of quenchings for H x > 0. We will see that the general theorems constrain this number for 
any value of J, and then show how this same conclusion comes about from the instanton 
analysis. 

We conclude the paper in Sec. with a summary and discussion. 

II. SU(2) INSTANTON FORMALISM REVISITED 

The instanton method is an efficient way of calculating tunnel splittings, both for parti- 
cles and for spin |26| . It is based on evaluating the path integral for a certain propagator 



in the steepest-descent approximation, and is designed to be asymptotically correct in the 
semiclassical limit (J — > 00 or h — > 0). Instantons are classical paths that run between 
degenerate classical minima of the energy. By "classical" we mean that the path obeys 
the principle of least action, and, a fortiori, satisfies energy conservation. However, a path 
along which energy is conserved can not run between two minima and still have real coor- 
dinates and momenta. Hence, one must enlarge the notion of a classical path, and allow 
the coordinates and/or momenta to become complex. One is naturally led in this way to 
reexamine the least action principle. In the case of spin, there is yet another reason for this 
reexamination, as we shall see. 

A. Spin coherent states and their parametrization 

More specifically, let the degenerate classical energy minima be along the directions rij 
and hf. We seek the propagator 

K fi = (h f \exp[-HT\\n i ) (2.1) 
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in the limit T — > oo. Here, |ni,/) are spin coherent states (defined below). This propagator 
is given by the path integral 

K fi = J V[h(t)}e- S ^\ (2.2) 

where n(0) = rij, n(T) = hf, and S is the (Euclidean) action. When we evaluate this 
integral by steepest descents, least action paths emerge naturally. 

We define the spin coherent state |n) as the state with maximal spin projection along the 
direction n. In other words, n is an eigenstate of J • n with eigenvalue J: 

J-n|n) = J|n). (2.3) 

The most common way of explicitly writing |n) and 5*[n(r)] is in terms of the spherical polar 
coordinates 9 and of the direction n. For formal purposes, another representation is more 
convenient. Let z be a complex number, related to (9, </>) by the stereographic map 

Q 

z = tan-e^. (2.4) 



Then, up to normalization and phase, |n) is identical to the state f27j 



\z) = e zJ ~\J, J). (2.5) 

The advantage of this representation is that matrix elements of various operators have nice 
analyticity properties. For example, for Hamiltonians polynomial in the components of J, 

= W (2 ' 6) 

is holomorphic in z and antiholomorphic in z', where z is the formal complex conjugate of 
z. In equations, this means that 

±H&,z) = ±H&,z)=0. (2.7) 

Explicit model calculations, however, are often easier in 9 and 4>. 

Since we will be considering paths on the complexified unit sphere, it is useful to under- 
stand what this means in terms of the two coordinate systems. Let u, v, and w be Cartesian 
coordinates in three-dimensional space. The real unit sphere is the surface specified by 

u 2 + v 2 + w 2 = l. (2.8) 
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The complexified unit sphere is obtained by allowing u, v, and w to become complex. The 
real and imaginary parts of Eq. ( p.8| ) provide two constraints among six variables (the real 
and imaginary parts of u, v, and w), leaving us with a four-dimensional manifold. We can 
also see this in the terms of polar and stereographic variables. Consider the former first. We 
relate them to Cartesian coordinates in the usual way 

(u,v,w) = (cos 9, sin 9 cos 0, sin 9 sin 0) . (2.9) 

Equation ( |2.8| ) is automatically satisfied. This continues to be true if we allow 9 and <fi to be 
complex. Once again, we conclude that four real quantities are required to specify a point 
on this manifold. Let us consider stereographic variables next. Since 

9 9 

z = tan - e** z = tan -e~^ } (2. 10) 

we see that if 9 and <p are complex, z will not be the same as z*, the true complex conjugate 
of z. To specify a point on the complex unit sphere, both z and z are needed, i.e., four real 
parameters are needed. Conversely, a point with stereographic coordinates (z, z) lies on the 
real unit sphere (which is a submanifold of the complex unit sphere) if z and z are complex 
conjugates. Such points may be given by specifying z alone or z alone. Thus, corresponding 
to the directions and hf, which are real, we may speak of "the points" Z( and Zf, or "the 
points" 2i and Zf. 



B. Action functional for spin 

Let us now consider the action. To specify a path n(t) on the unit sphere stereographically, 
one must give both z(t) and z(t). The action (or more precisely, the Hamilton principal 
function) for a path obeying the boundary conditions 



z(0) = Zi, z(T) = z f , 



(2.11) 



is given by p8|, p9| 



S(z{t),z{t)) = S K + S D + S_ 



(2.12) 



where 



S 



K 



J- 



zz — zz 



l + zz 



dt 



(2.13) 
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S D = f T H(z,z)dt (2.14) 



o 



Sb = —Jim 



{l + z{0) Zi ){l + z f z{T)) 



(2.15) 



{l + Z*Zi)(l + Z* f Z f ) 

We refer to these terms as the kinetic, dynamical, and boundary terms, respectively. In 
Sb, z*j denotes the true complex conjugate of Zij. This term depends explicitly on the 
boundary values of the path, and is needed to avoid the overdetermination problem [IB] . 
We state why this is so without proof. If we vary the path z(t), z(t), including the endpoints, 
and require the resulting variation 5S to vanish, we discover (i) the Euler-Lagrange (EL) 
equations 

^ (l + zz) 2 dH . (l + zz) 2 dH 

(2.16) 



2,/ dz ' " 2J dz ' 

and (ii) that SS has no terms proportional to Sz(0) or Sz(T). This means that the action 
evaluated for the classical path z c \(t), z c i(t) [which obeys Eqs. ( |2.11 ) and ( E-lGj )] is not a 
function of z^ = z(0) and Zf = z(T): 

-?-S d = -^-S cl = 0. (2.17) 

OZi OZf 

Equivalently, one can say that 

S[z cl (t),z d (t)] = S c \z f} z l} T), (2.18) 

where we show the variables on which S cl depends explicitly. Third, one obtains the 
Hamilton- Jacobi equations 

27 m ds« -z(o) 

dz f ~ l + z f z{TY dz t ~ J l + z(0) Zl - [ ' 

Lastly, the "energy" is conserved along the classical path: 



dt 

This follows from Eq. ( p. 16 ). 



^-H(z d (t),z cl (t)) = 0. (2.20) 



What would have happened if we had omitted the boundary term Sb in Eq. (|2.12| )? 
Since a general variation SS would include terms proportional to all four quantities 5zi, 5z~i, 
Szf, and 8zf, the classical path would have to be specified via four boundary conditions 
(zi, z~i, Zf, Zf), whereas the EL conditions would still form a system of differential equations 
of order 2. The system would be overdetermined with no solution in general. (The same 
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problem can be seen in spherical polar coordinates. One must give 0$ and 0j to specify n^, 
and Of and 0/ for hf, i.e., four boundary conditions in all.) Second, even with only two 
boundary conditions z(Q) = Zi, z(T) = z~f, for general Zi and Zf there is no real solution 
to the EL equations. That is to say, z c \(t) is not equal to the true complex conjugate of 
z c \{t). A solution will in general exist only if we allow z c \(t) and z c \(t) to be independent 
complex functions. [In spherical polar coordinate language, we must allow both 9(t) and 
(f>(t) to become complex.] We thus see that if we want to consider classical dynamics for 
spin from the viewpoint of a least action principle, we must allow the dynamics to live on 
the complexified unit sphere from the very outset. 



We have couched the above discussion in language following Faddeev [15]. It is also useful 



to discuss Klauder's resolution of these issues ||16|| . On the face of it, his treatment is slightly 
different, but turns out to be equivalent to Faddeev's in actual applications. Klauder does 
not include an explicit boundary term in S, but argues that since the continuum path integral 
is a formal construct with meaning only as a limit of its discrete version, one may add a 
term to the integrand for Sk that is quadratic in the velocities 9 and <f>, and which has an 
infinitesimal coefficient e. The EL equations are then a fourth order system, and a classical 
solution always exists. However, it has the following structure. It evolves from (9i,<pi) to a 
point (9i,4>i) (note the different use of the overbar) in a boundary layer of duration 0(e), 
evolves for a time T — 0(e) along a path 9(t), <f)(t) according to 

dH . d<b dH 



iJ sin 9— 

dt 



i J sin 9 — 
dt 
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(2.21) 



and then evolves from a point (9f, <f)f) to (9f, <pf) in another boundary layer of duration 0(e) 
near t = T. [Note that Eq. (|2.21| ) is equivalent to Eq. ([2.16D .] The extra kinetic term in 
the action gives a non-zero contribution only from integration over the boundary layers, but 
this contribution is explicitly independent of e as e — > 0. The net classical action is then 
given by 

r 

cos |0j cos ~9f 



S' 



cl 



iJ(l — cos 
+2 J In 



+ H(9, 



dt 



cos \9i cos ~9f 



[2.22) 
[2.23) 



The last term is the boundary layer contribution, and is equivalent to the explicit term Sb 
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in Eq. ( |2.12| ) ||30|| . The boundary values (#«,0i) and {0f,4>f) are constrained by 

tanf^e* = tan|^e^, (2.24) 
tan^/e-^/ = tan^e"^, (2.25) 

[note the similarity to z(0) = Zi, z(T) = zf\. Once again a solution can in general be found 
only if 9{t) and 4>(t) are complexified. 



C. Structure of answer for Fes 

Once the action functional is specified, the instanton recipe for calculating the tunnel 
splitting is as follows. Let there be a number of instantons, i.e., least action paths, labeled 
by k, and let the actions for these various paths be Sf. The tunneling amplitude is given 
by 

A = $>,e-< (2.26) 

k 

This can be understood as arising from the path integral fl2.2|) in the following way. In the 



integral, the dominant paths (in the sense of the method of steepest descents), are those for 
which the action is stationary. These are responsible for the exponential factor exp(— Sf). 
The prefactor results from integrating out the Gaussian fluctuations around the fcth 
instanton. The full expression for Kfi must include a sum over multi-instanton paths, and 
this leads to a result which is essentially exponential in A, in a way which is now well 



understood 25 



Up to now, all explicit instanton calculations of tunnel splittings of which we are aware 
have the property that the initial and final points of the instantons lie on the real unit 
sphere. In other words, z(0) = z*, z(T) = z}. In Klauder's language, {0ij,4>ij) = {Qiji < t>ij)i 
so there are no boundary layers in the path. (At intermediate times, of course, the path is 
complexified.) This is due to simplifying special features present in the models studied. As a 
result, there is no need to include the explicit boundary terms in S or S , and the practice in 
all papers on spin instantons (including those written by one of us, AG), has been to forget 
about them altogether. In the present problem, however, we find that although instantons 
with end points on the real sphere still exist, one also has instantons for which this is not 
so. We refer to these as boundary jump instantons (BJI). It is essential to include the latter 
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FIG. 2: Schematic of instanton paths in Fes- We show a portion of real phase space (the real 
unit sphere) as a two-dimensional surface, that is embedded in complex phase space. The latter 
is four-dimensional, but we can only represent it as three-dimensional in this perspective drawing. 
The initial and final points i and / lie in real phase space. Paths 1 and 2 start at these points, 
have no jumps, interfere with each other, and evolve smoothly as the C term in Eq. ( |1.4| ) is turned 
on. Except at the end points these paths lie in complex phase space. Paths 3 and 4 possess jumps 
at the end points, do not interfere, and are obtained only when 0. These paths lie entirely in 
complex phase space. 

in order to understand why the number of diabolical points on the H x axis is reduced from 
10 to 4. 

More specifically, we have four instantons, labeled 1, 2, 3 and 4 (see Fig. |2|). The first 
two have no boundary jumps, and are the ones that interfere as in Ref. H. The third and 
fourth are boundary jump instantons. Hence, 

A = X> fc e"< (2.27) 

k=l 

By proper choice of gauge, one can ensure that all Di are real, Sf and Sf are real, and that 
Sf = (Si 1 )*, D<i = D\. More generally, the real parts of Sf and Sf must be equal to each 
other by symmetry, and the imaginary parts must be related by 

Sf - Sf = 2zJ$, (2.28) 
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FIG. 3: Real parts of the actions versus the magnetic field for the four instantons in the Fe$ 
problem. Also shown is the energy splitting between the ground level pair. For h > ho, instanton 
number 3 is the dominant one, and since it has no interfering partner, the splitting rises with h 
instead of oscillating. 

where $ is half the area on the complexified unit sphere enclosed by the closed loop formed 
by taking instanton 1 from fij to n/, and instanton 2 back to hi. Therefore, we can write 



All the quantities in this equation depend on the field H x . We discover that Sf < Sf for 
all fields. At low fields, Re5f < Sf, while at high fields, the inequality is reversed. The two 
actions are equal at ho = 0.25 (see Fig. fj). The prefactors Di set the dimensional scale for 
the tunneling, and are generally equal to classical small oscillation frequencies in order of 
magnitude. Thus we do not expect the Di to be very different, and the relative importance 
of the different instantons is determined mainly by the Si. Hence, ignoring a very small 
region in the immediate neighborhood of ho, we can write 



In particular, there will be no quenchings in A for h x > ho. For h x < ho, quenching will 
occur when $ = (2n + l)7r/2J, where n is an integer. 



A = 2D x e 



-ReSf 1 



COS ( J$) + D 3 e~^ + D 4 e 



(2.29) 
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III. ANALYSIS FOR FE 8 



A. The interfering instantons 

We now turn to finding the instantons explicitly. For this purpose, it is better to use 
a coordinate system having z as the hard axis and x as the easy axis. Introducing spher- 
ical polar coordinates 9 and in the standard way, the energy (expectation value of the 
Hamiltonian) in this frame is given by 

E(6,4>) = Asin 2 flsin 2 + cos 2 #-2/icos# 

-2A 2 (cos 4 9 + sin 4 9 sin 4 0-6 sin 2 9 cos 2 9 sin 2 0). (3.1) 

The first step is to find the minimum of this energy. Setting dE/d9 and dE/dcf) to zero, 
we obtain 



cos sin ( 



24A 2 sin 2 #cos 2 #-8A 2 sin 4 #sin 2 + 2Asin 2 # =0 (3.2) 



and 



sin0[2Acos0sin 2 - 2 cos 9 + 2h + 8A 2 cos 3 9 

- 8A 2 (sin 2 9 sin 2 + 6 sin 2 9-1) cos 9 sin 2 0] = 0. (3.3) 

When these conditions are examined carefully, it is found that the minima occur when 
= 0, 7r and the expression in the square brackets in the second equation is zero, i.e., at 
= 0, Ti and at 9 = 9 where 9 obeys 

cos#o-^-4A 2 cos 3 fl = 0, (3.4) 

The minimum energy is 

Emm = cos 2 6*o — 2h cos 6*o — 2A 2 cos #0, (3.5) 

Since ( = 4A 2 /i 2 1 for all h, one can solve Eq. (|3.4| ) perturbatively to get cos^o = h+4X 2 h 3 
to first order in (. In the same approximation, E m m — —h 2 — 2A 2 /i 4 . 

Next, let us find the instantons. The trajectory, i.e., the path in phase space without 
regard to the time dependence, can be found by exploiting energy conservation. With the 
abbreviations 

u = cos9, s = sin0, (3.6) 
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FIG. 4: Plot of g(u, s) versus u for various s. The points to note are that (i) except when s = 0, 
g has no zeros in the interval [—1, 1], (ii) g always has two zeros outside this interval. For s = 0, g 
has a double zero at some u £ [—1,1]. The inset shows an enlarged view of g in the same interval. 

the condition E(6, 0) = E min can be written as 

g(u, s) = -^Z^u 4 + R(s)u 2 - 2hu + W{s) = 0, (3.7) 

where 

Z(s) = 4A 2 (l + 6s 2 + s 4 ), (3.8) 
R(s) = l-As 2 + 12A 2 s 2 + 4A 2 s 4 , (3.9) 
W(s) = g + h 2 + As 2 - 2A 2 s 4 , (3.10) 

with 

g = -(A + h 2 ) - E min « 2X 2 h 4 + 0{h Q ). (3.11) 

Equation ( p.7| ) may be solved as a quartic equation for u{s). Let us first consider only 
real values of <f), and thus of s. It is useful to sketch g(u, s) as a function of u for real u, 
and fixed s, remembering that \s\ < 1 (see Fig. f|). Consider first the interval — 1 < u < 1, 
corresponding to points (9, (p) on the real unit sphere. Since g = E — E min , it must be 
nonnegative in this interval. In fact, for s ^ 0, it must be strictly positive, and for s = 0, 
it vanishes only at one point, u = cos8 , where it has a double zero. Secondly, g — > — oo as 
u — > ±oo. It follows that for \s\ < 1, g(u,s) always has exactly two real roots u(s), with 
\u(s)\ > 1. These roots can not be the instantons that tend to the true energy minima as 
t — > ±oo. Those roots are the complex conjugate pair with both real and imaginary parts. 
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For the complex roots let u = A + iB with A and B being real. From the imaginary part 
of Eq. ( |3.7| ) we obtain the equation 

- 2Z(A 3 B - AB 3 ) + 2RAB - 2hB = 0, (3.12) 

while the real part yields 

- \Z(A 4 - 6A 2 B 2 + B A ) + R(A 2 - B 2 ) - 2hA + W = 0. (3.13) 
Since we are not interested in solutions with either A = or B = 0, Eq. ( |3.12| ) implies that 

B 2 = A 2 - (R/Z) + (h/AZ). (3.14) 
Substituting this into Eq. (|3.13| ), we obtain an equation for A alone: 

AZ 2 A 6 - ARZA 4 + (R 2 + 2WZ)A 2 - h 2 = 0, (3.15) 

which is a cubic in A 2 . 

When A 2 = 0, Z(s) = 0, and Eq. has the solution A = h/R = h/(l - A sin 2 0). We 



seek that solution of the cubic which tends to this solution as A 2 — > 0. We can obtain this 
approximately if we assume that A = O(h). The the terms ZA 4 and Z 2 A 6 are of order ( 
and ( 2 respectively, relative to the remaining two terms in Eq. ( |3.15| ). If we drop the former 
two terms, the remaining equation is trivially solved to obtain 

A= (R 2 + 2\vzy/2> (3 - 16) 

which is in fact of O(h). Thus, our assumption is self-consistently verified, and the solution 
has the correct behavior as A 2 — >• 0. 
Next, we note that 

R 2 + 2W Z = Pq + P' 2 sin 2 (f) + P[ sin 4 0, (3.17) 

where we have unabbreviated s, and 

p> = i + 8A 2 (/i 2 + 2o)~(l + C) 2 , (3-18) 

P' 2 = -2A + 24A 2 + 8AA 2 + 12( + 48A 2 # , (3. 19) 

ss -2A + 24A 2 + 8A 2 A + 12C + 6C 2 , (3.20) 

Pi = A 2 + 8A 2 + 24AA 2 + 128A 2 + 2( + 8A 2 # - (3.21) 

w A 2 + 8A 2 + 24AA 2 128 + A 2 + 2C + C 2 - (3-22) 
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FIG. 5: Plot of Re cos 9 versus 4> for h = 0.2 and the Fes parameters. The solid line is from exact 



numerical solution of g(u, s) = 0, and the dashed line is the approximation ( 3.23 ). 



All three coefficients depend on h only through the combination ( = 4A 2 /i 2 , which is very 
small. If we neglect this weak dependence, we get A pa Aq, where 

h 

A ° = (1 + P 2 sin 2 + P 4 sin 4 0) 1/2 ' (3 ' 23) 

with 

P 2 = -2A + 24A 2 + 8A 2 A, (3.24) 

P 4 = A 2 + 8A 2 + 24AA 2 + 128A 2 . (3.25) 

For completeness, we mention that it is possible to systematically obtain corrections to 
A in powers of (, in the form A Q + (Ai + ■ • -. We do not carry out this exercise here. To 
see how good Eq. ( |3.23| ) is, we plot the real part of cos 9 in Fig. |5|. The dashed line is our 



approximation, and the solid line is obtained by numerically solving E — E m \ n = with the 
numerically exact value of E m \ n . In this plot, h is taken to be 0.2. It can be seen that the 
two curves agree very well. 

With Eq. ( 3.23|) , we can now calculate the imaginary part of the tunneling action, and 



the Berry phase 2J$ of Eq. (|2.28|) . We have 

$ = r(l-A)d(j) 
Jo 

pa [7r-hI(X,X 2 )}, (3.26) 

where 



jo (1 + P 2 sin <p + P 4 sm <p) 1 i z 
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FIG. 6: Tunnel splitting (in Kelvin) between ground level pair for the model (L4), as computed 
by numerical diagonalization of the Hamiltonian (solid line) , in the instanton approximation keep- 
ing only instantons 1 and 2 (dot-dashed line), and in the instanton approximation with all four 
instantons (dotted line). We take all prefactors Di to be equal and independent of h, and adjust 
the common value so as to obtain the correct answer for A at h = 0. 

Various formulas for this integral are given in the Appendix. 

The quenching condition J$ = (n + \)tt gives the diabolical points as 

(2J-2n-l)7T 



2JI 



[3.28) 



where n is an integer. For Feg, / = 3.88 as stated in Sec. I. The observed diabolical 
points agree extremely well with Eq. ( |3.28| ). In particular, Eq. ( |3.28| ) gives the period 
AH = 0.409 T. The experimental period is 0.41 T. 

We can now see why the diabolical points in Fes are so regularly spaced. This is because A 
is linear in h to very good approximation. Corrections to Eq. ( |3.28|) can be found by including 
the correction A\. This may be important for the new system in which A oscillations are 
indicated 0. 

At this point, we could use Eq. (|3.14j) and our approximation A « Aq + C^4i to find B((p), 
and thus ReS^ 1 . Since the resulting analysis is not completely analytical, we forgo it, and 
instead, solve for the instanton trajectory cos (9(0) and evaluate the integral for the action Si 
numerically. The instantons are shown in Fig. 3 of Ref. [D], and the resulting approximation 
for A is shown in Fig. |6|. We do not know the prefactor Di, but it is clear that it is a very 
good approximation to take it to be independent of H . It is also clear, however, that the 
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interfering instantons can not account for the behavior of A for h > ho- For that, we must 
turn to the boundary jump instantons. 



B. The Boundary- Jump (Noninterfering) Instantons 

The instantons we have found above have the properties (a) Zi = z*, Zf = (%)*, and (b) 
H(zi,Zi) = H(zf,Zf) = E min . As discussed earlier, however, only Z; L and Zf are fixed, and 
one need not have Zj = z* or Zf = (%)*. Let us suppose that we start the evolution of the 
Euler-Lagrange equations ( [2.161 ) from a point (z i7 z^, with Zi ^ z*. Whatever this point is, it 



follows from Eq. ( |2.20| ) that the energy H(z, z) must be conserved along the trajectory, and 
must be equal to E — H(zi, z{). Conversely, if the value of E is known, we can determine 
the possible values of z\ by solving the equation H(z, zi) = E. For the non-boundary- jump 
instantons, E = E m i n . The issue is what value of E we should use for the instantons with 
jumps. 

The answer is that we must take E = E m[n for the jump instantons too. The easiest way 
to see this is in Klauder's formalism. His extra kinetic term is 

s> K = Ae f rrr- \2 mt - ( 3 - 29 ) 
Jo (l + zzy 

It is easy to write down the new Euler-Lagrange equations, and see that energy is once again 
conserved. Suppose there is a boundary layer in the solution of the EL equations around 
t = to- Since, as e — > 0, the boundary layer turns into a jump, and the extra kinetic energy 
vanishes for times t = t + or t — , the energy before and after the jump must be the same. 

The problem of finding the boundary-jump instantons (indeed, all instantons) can there- 
fore be posed as follows. Let the "classical" Hamiltonian H(z',z) have minima at (z*,Zi), 
(zf, Zf), and let its value at these points be E min . Then we find all possible values of Zj by 
solving the equation 

H(z, Zi) = E min . (3.30) 
This equation has a double root at z = z*, since 



0. (3.31) 



It may, however, have additional roots at z^ ^ z*. These will be the initial points of 
boundary jump instantons. (Final points Zf are obtained in exactly the same way by solving 
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H(zf,z) = -Emm-) We then find all possible instantons z(z) as solutions to the equation 
H(z, z) = -E m in, and identify which solution connects on to which end point. The time 
dependence is not needed to compute the action. For, the kinetic term can be written as 



S 



K 



-J 



1 + zz(z) 



dz 
dz 



ziz 



dz 



r4> } 
%J / 1 



cos 



the dynamical term equals 



S 



D 



E 



(3.32) 
(3.33) 

(3.34) 



for all instantons, and the boundary term Sb depends only on the boundary values, vanishing 
for the instantons without jumps. 

For our Fes Hamiltonian, the argument following Eqs. (|3.7| )- (|3.11| ) shows that we have 
already found all the instantons without jumps. To find the jump instantons, let us first 
note that (0j = 9q, fa = 0), and (Of = 6q, 0/ = ir). The constraints on Z{ and Zf therefore 



reduce to 



tan~0je* = tan|0 o > 
tan-flye - **' = -tan§0 o . 



(3.35) 
(3.36) 



Because of the symmetry of the problem, solving either of these equations will be enough 
for deducing the solution of the other. Let us solve for the initial values. First, Eq. ( |3.35| ) 
can be solved for sin fa to yield 

• 2 7 ( COS0j - COS O V ,a QV x 

sm fa = -[ . - , . (3.37) 

y sm &i sm u J 

Substituting this formula and Eqs. flOp and ( p.5|) into the energy conservation condition 
E(0i, fa) = E min yields 







H 

2A 2 



cos 9j — cos Go 



sin 
cos 6f + 



+ (cos 2 6i — cos 2 6* ) — 2/i(cos 9i — cos 9 



'cos6>i - cos# n \ 4 _ 2 x / cos 9i - cos \ 2 4/1 
O COS 9i I ; — ; — COS 9q 



.(3.38) 



sin^o / V sin^o 

We now eliminate h from this equation using Eq. ( |3.4| ). After some straightforward but 
lengthy algebra, we obtain 

A 



(cos^j — cos^o) 2 



1 - 



sin 2 9r] 



2X 2 b(9 i7 9 ) 



0. 



(3.39) 
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where 



b(9i, ) = cos 2 9i + 2 cosfc cos9 + 3 cos 2 9 + V ; A , , + G^TT- (3.40) 



(cos 9i — cos 9 ) 2 gCos 2 0j 
sin 4 sin 2 9 



There are four solutions to Eq. ( |3.39 ). The first two are the non-jump solutions, cos Qi = 
cos^o- The other two, which are the jump solutions obey 

l--^--2X 2 b(9 t ,9 )=0. (3.41) 
sin 9q 

This is a quadratic equation for cos 6^. For h < 1 and small A2 it is easy to check that in 
both solutions, cos Oi is real and greater than unity. Hence, we may write 9{ = ivq, where vq 
is real. Equation ( 3.35 ) then shows that 0, may be taken in the form ^tt — ifi , with fi real. 

It is easy to see that Eqs. (|3.38|) -( 3.41 ) continue to hold if 9i is replaced by Of. Thus the 
possible values for Of are the same as those for 0^. Numbering the instantons in question 3 and 
4, we have either 0^ = 0f \ or 0^ = 0f \ Symmetry suggests (and explicit numerics verifies) 
that the latter possibility is the correct one. If we then divide Eq. (|3.36| ) by Eq. ( |3.35|) , we 
see that 0/ = tt — 0j. Thus, the end points are related by the symmetry of reflection in the 
hard-medium plane, i.e., (J e , J m , Jh) — > (—J e , Jm, Jh), where the suffixes e, m, and h stand 
for easy, medium, and hard. 

If we parametrize 0(f) and 0(t) as iv(t) and (tt/2) + ifi(t), then it easy to verify from 
the equations of motion that fi(t) and v(t) stay real at all t. With this parametrization, the 
kinetic term in the action for instantons 3 and 4 may be written as 



/"MO 

S K = J [coshu(fi) - l]dfj,, (3.42) 

J— IXo 

which is real. The boundary contribution is 

cosh(z/ /2) 



S R = 4 J In 



(3.43) 



cos(fl /2) 

which is also real. Hence the total action for the jump instantons is real. 

The explicit calculation of the actions must be done numerically. We solve for the tra- 
jectories in the form using energy conservation, making sure that the end points are 
correct. All these calculations are done as a function of h. The results for have already 
been shown in Fig. |3|. We can also calculate the splitting using Eq. ( |2.30| ), taking D 3 = Di, 
and fixing D\ as before. The result is shown in Fig. |5|. As can be seen, the agreement with 
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FIG. 7: Same as Fig. ||, but for C = 1.2 x 1CP 5 K. In addition to the splitting A obtained by 
numerical diagonalization of the Hamiltonian matrix (solid line), we show the answer given by 
Eq. ( |2.29D (dashed line) with all prefactors chosen to be equal, and adjusted so as to agree with 
the numerically computed splitting at h = 0. 



the exact diagonalization is rather good. The t-dependence of J for instanton number 3 is 



shown in Fig. 2 of Ref. [17]. The jumps are evident in this figure. 

We have also carried out this exercise for C = 1.2 x 1CT 5 K. The results are shown in 
Fig. [7]. The general quality of the instanton approach is again very good, but it is weaker 
near the point where Sf = Re Sf. The obvious reason is that we have not considered the 
variation with h of the prefactors, especially D 3 . If we decrease C by yet another factor of 
4, Eq. ( |2.29| ) with a single, /i-independent prefactor D 1 = D 3 , gives six quenching points 
instead of eight as found numerically. These considerations show that the prefactors are not 
always unimportant; at present, however, we only know how to find them for the non-jump 
instantons [pTj], and their calculation for jump instantons is an open problem. 



C. Is the fourth-order anisotropy a singular perturbation? 

In the previous subsection we showed that jump instantons exist for any A2 7^ 0. The 
character of the least action or "classical" paths for the problem is qualitatively altered by 
the fourth order term in Eq. ( |1.4| ), and from the point of view of the semiclassical analysis, 
therefore, this term is a singular perturbation. Viewed as a quantum mechanical operator, 
however, it is clear that the term is nonsingular; an infinitesimal nonzero value of A2 can not 
change the qualitative nature of the energy spectrum. We therefore refer to the perturbation 
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as quasisingular. 

It follows from this consideration that a formal analysis of the instanton approach in the 
A2 — > limit will necessarily be rather delicate. We describe briefly some analysis that shows 
why this is so. 

We focus on jump instanton number 3, as this is clearly the most important new contrib- 
utor. For this instanton, we discover that 

^0 = ~ln(^)+O(A 2 )V2, (344) 



where uq = cos#o 5 which is given by Eq. ( |3.4|) . Secondly, to leading order in A2, 

1/2 

(3.45) 



cosh v 



V^2 



1 — A cosh 2 n 



1 + 6 cosh 2 jj, + cosh 4 jj, 



Because of this A 2 ^ 2 dependence, the leading term in the action Sf is also of order A 2 
The next term is proportional to JlnA 2 . While the latter is well behaved as A 2 — > 
(remember that S cl must be exponentiated), the former is not. It is clear that there must be 
a cancellation due to a corresponding term in D3, but we have not attempted to find this. 
This analysis shows that the success of our assumptions about the scale and /i-independence 
of the prefactors is somewhat fortuitous. However, this assumption is the natural one in a 
semiclassical approach. 

IV. NUMBER OF QUENCHING POINTS 
A. Berry phase argument 

That the tunnel splitting between two states vanishes is just another way of saying that 
the states are degenerate. The magnetic fields at which such degeneracy occurs form a set 
of isolated points in the magnetic field space (H x , H y , H z ), that are said to be diabolical, 
following Berry and Wilkinson [ 32f . Since these points are singularities of the energy surface, 



there are strong constraints on their creation or destruction as a perturbation is continuously 
varied. The number of diabolical points on the H x = axis is known when C = 0; it is 
interesting to inquire how many may remain when C is turned on, and to pursue this inquiry 
for general J, not just J = 10, and also consider fields in the xz plane f33|. (Simple physical 



arguments show that there can be no degeneracy if there is a component of the field along 

y-) 
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FIG. 8: Trajectories of diabolical points under the influence of the fourth order perturbation 
— C{J\ + </£) in Eq. ( |1.4| ). The numbers and ±7r are the Berry phases associated with the 
adjacent contours. 

Let us first take C — 0. With H y ^ 0, in the standard representation of the spin 



operators, the matrix of the Hamiltonian is real. By a general theorem [33], the codimension 
of a degenerate eigenvalue of a real Hermitean matrix is 2. Hence, the degeneracies must 
occur at isolated points in the H x -H z plane. When C is turned on, each one of these points 
must turn into a line in the three-dimensional (H x , H Z ,C) space. The only two kinds of 
behaviors that are permitted by the theorem are shown in Fig. || The first kind, marked 
'a', is a diabolical point that continues on for ever. The second kind, marked 'b', shows that 
two distinct diabolical points in the C = plane actually lie on the same diabolical line in 
the three-dimensional space. (Of course, a similar turnaround could connect the two lines 
marked a and b for some C < 0, or the line b could form a closed loop.) 

Is it possible for one of the diabolical lines to terminate abruptly? An argument based on 
Berry's phase |35] shows that the answer is no. Let H now denote the two-dimensional vector 
(H X ,H Z ). Let two states \ipi(C, H)) and [^(C, H)) be degenerate at H = H = (H x0 ,H z0 ) 
for some value of C, and let Y be a small closed contour in the H x -H z plane around the 
point Ho- Berry's phase, given by 

7 (r) =i jf (V>i(C,H)|Vh^(C,H)) s (4.1) 

equals ±7r if Y encloses a diabolical point. Otherwise, 7 = 0. Since a small change in C or H 
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gives rise to a nonsingular perturbation of the Hamiltonian, the state ipi(C, H) is a smooth 
function of C and H. Hence the integrand of Eq. ( |4.1| ) can not change discontinuously under 
a continuous change of C, and the integral must not change at all. Thus the contour T must 
continue to encircle a degeneracy at small non zero C if it did so at C = 0. 

From this point of view, the behavior 'b' in Fig. ^ can arise only if 7 has opposite values 
for the two diabolical points at C = 0. The Berry phase for a contour r 2 encircling both 
points is then 0, and it is possible that for larger values of C, the contour can be shrunk to 
a point without running into any singularity. Plainly, this can happen only if the two points 
annihilate each other at some C. We can think of simply slipping the contour r 2 off the 
diabolical line by lifting it above the hairpin bend in the figure. 

It follows that diabolical points can only disappear in pairs. For the problem of interest 
to us, degeneracy between ground levels, the points are constrained to occur when H z = 
(consider the behavior of the Hamiltonian under a 180° rotation about x). Our analysis 
shows that with increasing C, the points at larger H x are removed first. Since the quenching 
points for any value of J are located as given by Eq. ( |1.2| ) when C = 0, it follows that when 
C^O, the number N q of such points for H x > must depend on J as follows: 

2J N q 

4n even 

4n + 1 even (4.2) 

4n + 2 odd 

An + 3 odd 

The same number must occur for H x < 0. And, if J is half-integral, there must be a 
quenching point at H x = 0, consistent with Kramer's theorem. For Fe§, J = 10 and N q = 4, 
which is consistent with these general rules. 

B. Instanton based argument 

Let us now see how our instanton analysis yields the same conclusion. For this, let us 
ignore the fourth instanton, and write the amplitudes due to the remaining instantons as 

H 

A 1+2 = Dxie-^ + e- 3 ^), (4.3) 
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A 3 = D 3 e- s Z . (4.4) 

The keystone of our argument is the relative sign of the two amplitudes A 1+2 and A 3 , in 
particular, the fact that 

sgn(A 1+2 ) = sgn(A 3 ) for large h. (4.5) 

By "large h" , we mean that h is just less than the field strength at which the two classical 
minima in the energy merge into one. We shall not try and prove Eq. (|4.5|) with mathematical 



rigor | 57| . Rather, we argue that it is physically plausible, for at such large fields, only spin 
orientations in a very small angular range are important, and one can use the Holstein- 
Primakoff, Villain, or any of a number of similar mappings to approximate the spin operators 
in terms of Q and P, position and momentum operators for a particle in one dimension, The 
problem is thereby mapped on to a particle in a double well in one dimension, for which 
the splitting never vanishes. Were Eq. ( |4.5| ) not true, it is conceivable that we could make 
A 1+2 + A 3 vanish by varying the relative value of C and (k\ — fc 2 ). 

The second point is that the sign of A 3 should not change with h, since instanton 3 acts 
alone and has a real action. We may therefore take _D 3 > 0, so that A 3 > for all h. 

The third and last point is that the sign of Ai + 2 at h = is now fixed by the requirement 
that this amplitude vanish the correct number of times between and large h. Readers can 
verify that a correct assignment is obtained by taking 

Ai = 2D ie ~ RcS ? cos(Jvr), {h = 0), (4.6) 

with Di > 0. 

The rest of the argument is simple. Continuing to ignore the fourth instanton, the 
condition A = can be rewritten as 

A 3 = -A 1+2 . (4.7) 

We now simply sketch both sides of Eq. (|4.7|) as a function of h, keeping in mind the three 
points made above. This is done in Fig. |9] for all four classes of J listed in Eq. (|4.2| ). In 
each case, it is obvious that the number of zeros, N q , is exactly as given in this equation. In 
particular, it is obvious that N q can change only in steps of two if the curve for A 3 is raised 
or lowered [[38 . 
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FIG. 9: Sketch of — A1+2 and A3 versus h for all four classes of J listed in Eq. ( |4.2| ). The 
quenching points are given by the intersections of these two curves. The key points to note in each 
case are the number of zeros of Ai + 2 and the number of quenching points. 

V. CONCLUSION 

We have shown that the instanton formalism for spin coherent state path integrals re- 
quires the inclusion of instantons with discontinuities at the end points as a general matter. 
Such instantons are essential to understanding the magnetic field dependence of the tun- 
nel splitting in Fes- We have shown that with certain plausible assumptions about the 
prefactors, the instanton approximation can be quantitatively accurate. However, proper 
calculation of tunneling prefactors for the instantons with jumps remains an open problem. 

Since jump instantons arise as a result of overspecification of boundary conditions, and 
since this overspecification is a necessary consequence of the coherent state formulation, 
it is clear that similar instantons will in general be present in all path integrals based on 
coherent states. For spin, such path integrals are unavoidable if one wishes to treat all spin 
orientations on an equal footing, and they are the only way of passing to the classical limit 
of the dynamics. This is not so for massive particles. Nevertheless, there may well be some 
problems that are better formulated in terms of coherent states, and then one will have to 
be alert to the presence of jump instantons. An explicit instance where this is so remains to 
be found. 

Acknowledgments 

One of us (AG) would like to acknowledge insightful discussions with Michael Stone on 
spin coherent state path integrals, and Jacques Villain on diabolical points. 



28 



* 



APPENDIX A: THE OSCILLATION PERIOD INTEGRAL 



In this appendix, we provide some formulas for the oscillation period integral fl3.27| ), 
which we reproduce here for convenience: 

J ( A ' ^ = f M-LP • JxP • 4^1/2 ^- ( A1 ) 

Jo (1 + P 2 sm <p + P 4 sm q>) ' 

For Fes, the constants P2 and P4 equal —1.200 and 0.7576 respectively. The integrand is 
real with these numbers. 

Let w be a complex number such that 

w + w * = -p 2) ww * = p 4 . (A2) 
Then, using formula 2.616.1 of Ref. | 55| , we have 



Jo [(l-wsin 2 0)(l-w*sin 2 0)] 1 /2 = ^]—^ 1-w ) ' ( ^ 

where PJ(/c) is the complete elliptic integral of the first kind. 

The form ( |3.27| ) shows that / is real, but this is not evident from Eq. ( |A3| ) . But, by using 
the canonical form for K, we obtain 

/■7r/2 rjnf- 

jo (1 — w cos^ x — w* sin xyi 1 

If we write w = a + ib, expand the integrand in powers of b, and integrate term by term, we 
obtain 



7T 



\/l — a 



b \\ , (-1)" / 6 y n r(2n+|)r(n+|) 



+••• + 



16\l-a/ 7T V 1 " / r(2n + 1) T(n + 1) 



VT=^~ V4'4'"' (1-a) 2 / 

where P is the hypergeometric function. This form again shows that I is real. With the 
Fes parameters, a = 0.5999, and b = 0.6307. The series ( |A5|) is not convergent for these 
values of a and b. For efficient numerical evaluation of / in such cases, one can either 
perform a numerical integration of Eq. ( |3.27| ) directly (which is what we did) , or reexpress 
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the hypergeometric function in Eq. (|A5|) in terms of other hypergeometric functions of the 
argument 

a-a) 2 = ( 2 + ^) 2 fA6l 
l-2a + (a 2 + 6 2 ) 4(1 + P 2 + P 4 ) 1 ; 

using formula 9.132.1 of Ref. |]39 |. Since the argument is now small compared to 1 (0.2870 
for Fes), the hypergeometric series involved are rapidly convergent. 
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